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Abstract 

At temperatures below the critical temperature of superfluid phase transition baryonic matter 



J3 i 

S^-f emits neutrinos by breaking and recombination of Cooper pairs formed in the condensate. The 

o 

strong interactions in the nuclear medium modify the weak interaction vertices and the associated 
neutrino loss rates. We study these modifications non-perturbatively by summing infinite series 
of particle-hole loops in the S- wave superfluid neutron matter. The pairing and particle-hole in- 
teractions in neutron matter are described in the framework of the BCS and Fermi-liquid theories 
derived from microscopic interactions. Consistent with the /-sum rule, the leading order contri- 
bution to the polarization tensor arises at 0(q 2 ) in the small momentum transfer, q, expansion. 
The associated neutrino emission rate is parametrically suppressed compared to its one-loop coun- 
terpart by a factor of the order of 5 x 10 -3 , the parameter being the baryon recoil in units of 
temperature. 
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I. INTRODUCTION 



Pair-correlated baryonic matter in compact stars emits neutrinos via the weak neutral 

nn 

current processes of pair-breaking and recombination 12J 

{NN} ^ N + N + u f + u f , N + N -> {NN} + v f + u f , (1) 

where {NN} refers to a Cooper pair, iV + A" to two quasiparticle excitations, Vf and \>f to 
the neutrino and anti-neutrino of flavor /. The process ([T]) is limited to the temperature 
domain T* < T < T c , where T c is the critical temperature of pairing phase transition and 
T* ~ 0.2 T c . At and above T c this reaction cannot occur, since momentum and energy cannot 
be conserved simultaneously in a process A" — > A" + Uf + Uf, i. e., an on-mass-shell fermion 
cannot produce bremsstrahlung (in the absence of external gauge fields). At asymptotically 
low temperatures, T < T*, the rate of the process (CQ) is exponentially small since the number 
of excitations out of the condensate is suppressed as exp(— A/T), where A(T) is the gap in 
the quasiparticle spectrum. 

Cooling simulations of neutron stars revealed the efficiency of the processes (CQ) in refrig- 
erating the baryonic interiors of a compact star from temperatures T < T c ~ 10 9 K down 
to temperatures of the order of 10 8 K The temperature domain above corre- 

sponds to the neutrino cooling era that spans the time domain 10 2 < t < 10 5 years. The 
predicted surface temperatures of neutron stars during this and the following photon cooling 
era (where the star loses its thermal energy by emission of photons from the surface) are 
sensitive to the neutrino emission rates within this time-domain. Remarkably, the process 
([1]) relates the cooling rate of a compact star to the microscopic physics of its interiors and is 
particularly sensitive to the density-temperature phase diagram of paired baryonic matter. 
Therefore, the measurements of the surface (photon) luminosities of neutron stars and their 

rave predictive power for analyzing the phase 



rave predici 



interpretation in terms of cooling simulations 
diagram and composition of baryonic matter 

The rate of the process ([I]) was computed independently (and within alternative methods) 
in Refs. in the case where the pairing interaction is in the 1 Sq partial wave, i. e., 

nucleons are paired in a spin-0, isospin-1 state, (the influence of electric charge carried 
by proton Cooper pairs and the case of pairs forming a spin 1-superfluid are discussed in 
Refs. Q). In propagator .anguage these rates correspond to the one-,oo P approximation 
to the polarization tensor of baryonic matter |2J, HQJ. It has been known for a long time 



that the gauge invariance requires summation of infinite series of loop diagrams (known 
as RPA resummation) in conventional superconductors 3, Q]- -F^f- 14 1 carried out RPA 
summations in neutron and quark matter within the context and kinematics of neutrino 
scattering in hot proto-neutron stars. Subsequently, Ref. [ijj] applied the gauge invariance 
and vertices derived by Nambu [12} to the neutrino emission processes in the vector channel 
and concluded that the rates are suppressed by by many order of magnitude. 

We shall arrive below at results which are at variance to those of Ref. [lj]] for the fol- 
lowing reasons. The effective vertices implemented in Ref. [lsj] were derived within a zero- 



temperature theory [l2j|; herein we shall derive the vertices at finite temperature, consistent 
with the polarization tensor of matter. This guarantees that the unitarity of the S'-matrix 
in th e q uantum mechanical process of the bremsstrahlung is preserved. Second, while in 
Ref. [15| the matrix elements are expanded in the parameter vf/c, where Vp is the baryon 
velocity, and the leading order contribution to the rates comes from the terms 0(vp/c) 4 , 
herein the polarization tensor is expanded in small momentum transfer, q; the leading order 
contribution to the rates arises at 0(q 2 ) consistent with the /-sum rule. As a consequence, 
we find a suppression of the neutrino bremsstrahlung rate which is by several orders of mag- 



nitude less than predicted in Ref. 



15] . Since the magnitude of the pairing gaps are not well 



known and the neutrino bremsstrahlung rates are rather sensitive to the value of the gap, 
the pairing braking process in S'-wave superfluids remains a potentially viable mechanism 
of cooling of neutron stars. 

We start by setting up a formalism for computing the vertex corrections to weak reaction 
rates that arise due to the strong interactions in baryonic matter. We consider the case of 
neutron matter at subnuclear densities, which we describe within the Landau Fermi-liquid 
theory derived from microscopic interactions. At these densities neutron matter is charac- 
terized by an isotropic order parameter arising form the interaction in the 1 S'o partial wave 
channel; we solve the corresponding problem of pairing in the framework of the Bardeen- 
Cooper-Schrieffer (BCS) theory. In the non-relativistic limit, the vector and axial vector 
weak vertices are associated with scalar and spinor perturbations. Thus, the problem of 
the weak vertex renormalization reduces to a study of the effective three-point vertices that 
sum-up particle-hole irreduceable ladders in superfluid matter in the scalar and spin chan- 
nels. In a wider context, such resummations describe of the low-frequency, long-wave-length 



collective modes of superfluid Fermi-liquids 13j, ll6|, Uj§. In particlar, such resummations, 
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known as quasiparticle random phase approximation (QRPA), have been widely used in 



nuclear structure calculations 



id- 



The paper is organized as follows. Section [TT] introduces the Green's functions formal- 
ism and computes, for the sake of illustration, the neutrino rates at one-loop. In Sec. II III 
the modifications of the weak interaction vertices are computed by summing particle-hole 
ladders in superfluid neutron matter. Section [TV] is devoted to the computation of the full 
polarization tensor that includes the vertex corrections derived in the preceding Section. 
Section [V] contains our conclusion. Some technical details are relegated to the Appendix. 

II. S- WAVE PAIR-CONDENSATION IN NEUTRON MATTER 



Below we shall describe the neutron pair-condensate at subnuclear densities within the 
framework of the Fermi-liquid theory, which assumes that the elementary degrees of freedom 
are quasiparticles with well defined momentum-energy relation and infinite life-time. The 
interactions between the quasiparticles are then described in terms of Landau parameters 
which depend on the momentum transfer (or scattering angle). Since the scattering angles 
involved are typically small, the momentum dependence of Landau parameters is further 
approximated by the leading and next-to-leading order terms in the expansion in Legendre 
polynomial with respect to the scattering angle. The problem of pairing in neutron matter 
will be solved below within the BCS approximation, where the anomalous self-energy (the 
gap function) is computed from the bare interaction while the normal self-energy is computed 
within the decoupling approximation which ignores the effects of pair-correlations on the 
single particle spectrum of quasiparticles. A number of factors such as the renormalization 



of the pairing interacti on Jl9l 
renormalization 
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24 



25 



26 



27 



281 ] and the wave-function 



25l . |29| . 130] . affect the absolute value of the gap. The role of these factors 



has not been settled yet, and we shall employ below the standard BCS approach which has 
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32 



led to convergent and verifiable results for neutron matter (see the reviews 

Since the baryonic component of stellar matter is in thermal equilibrium to a good ap- 
proximation, we shall adopt the Matsubara Green's functions for the description of the 
neutron condensate and for the evaluation of the polarization tensor. In the case of So 
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pairing these are defined as (e. g. Ref. [34J |. pg. 120) 

G aa '(p,r) = -8 aal (T T a pa (T)a pa ,(0)), (2) 
F aa/ (p,r) = 5 f? /(T r a_ pi (r)a pT (0)), (3) 
%*(P,t) = MT T at T (r)al pi (0)>, (4) 

where a stands for spin o =f, I and isospin r is the imaginary time, T T is the imaginary 

time ordering symbol, a) {r) and a pu {j) are the creation and destruction operators. In the 

momentum representation the propagators are given by 

/ u 2 v 2 \ 

G aa >{ip n ,p) = <W<W - h - — , (5) 

\ip n -e p ip n + e p J 

F aa >(ipn,P) = -icr v 5qq'u p v p ( - — — ] , (6) 

and Fl a ,(ip n , p) = F aa r(ip n , p), where p n = (2n+l)irT is the fermionic Matsubara frequency, 
a y is the ^/-component of the Pauli-matrix, u 2 = (1/2) (1 + £ p /e p ) and v 2 = 1 — u 2 are the 
Bogolyubov amplitudes and 



'e P + (?) 

is the quasiparticle spectrum, where £ p = p 2 /2m+ReS(p)— /z is the spectrum in the unpaired 
state with m and fi being the bare mass and chemical potential. Here £(j>) and A(p) are 
the normal and anomalous self-energies. For the later diagrammatic analysis we shall need 
the hole propagator which is given by 

Gi a ,(ip n ,p) = G aa <(-ip n ,-p) 

( u 2 v 2 \ 

= -6 m ,6tf(^2— + ^—). (8) 

\ip n + e P ipn-tpj 

The spin and isospin dependence of propagators for S'-wave spin-0 and isospin-1 pairing is: 
G aa '{ip n ,p) = 6 mr >6qq'G(ip n ,p) and F aa ,(ip n) p) = -ia y 5qq'F(ip n ,p). 

Since the quasiparticles are confined to the vicinity of the Fermi-surface we expand the 
normal self-energy around the Fermi momentum, pp, to obtain 

e(p) = ^(p-Pf) ~ H* , (9) 

m* 



m 
m 



771 

1 H <9 p Re£(»|p = p F 

Pf 



(10) 



where //* = — e{pp) + n — ReS(p^) is the effective chemical potential, m* is the effective 
mass. The dependence of the self-energies on the off-mass shell energy will be neglected, 




FIG. 1: Temperature dependence of the 1 5o pairing gap in neutron matter for Fermi wave vectors 
Pf = 0.4 (solid line), 0.8 (dashed line), 1.2 (dashed- dotted line) and 1.6 (dashed- double- dotted line). 

i. e. the wave-function renormalization which accounts for the next-to-leading term in the 
expansion around the Fermi-energy is set to unity. The BCS mean-field approximation to 
the anomalous self-energy can be written in terms of the attractive four-point vertex function 
T(p,p') in the form 

A(p) = -2 J ^fM) ^(p')f(u') , (11) 

where f{oS) = [1 + exp(/5u;)]~ 1 is the Fermi distribution function and (3 is the inverse tem- 
perature. For the problem of neutron matter at subnuclear densities we adopt the standard 
BCS approach and thus approximate the four-point vertex by the bare interaction in the 
1 S'o partial wave channel and the single particle spectrum by the on-shell spectrum given by 
Eq. (J9]). With these approximation Eq. ( ITTi) reduces to 



A(p) 



d 3 p' 



■V(p,p') 



A(p>) 



[1 " 2/(e(p'))] • 



2(2vr)3 V£(p0 2 + A(pO 

The gap equation is supplemented by the equation for the density of the system, 

d 4 p 



P 



(27T) 



-lmG(p)f(u), 



(12) 



(13) 



which determines the chemical potential in a self-consistent manner. Fig. [T] shows the tem- 
perature dependence of the 1 5'o pairing gap in neutron matter for several densities parame- 



TABLE I: Density dependence of the effective mass, the scalar and spin-spin interactions, the 
pairing gap and the critical temperature; the interactions are given in units of the density of states 
u(p F ). 
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[MeV] [MeV] 


0.4 


1.02 


-0.56 0.55 


1.54 


0.85 


0.6 


1.00 


-0.50 0.49 


2.60 


1.44 


0.8 


0.97 


-0.47 0.44 


3.15 


1.78 


1.0 


0.94 


-0.45 0.41 


3.09 


1.80 


1.2 


0.92 


-0.43 0.40 


2.44 


1.46 


1.4 


0.88 


-0.41 0.40 


1.41 


0.88 


1.6 


0.84 


-0.36 0.39 


0.57 


0.38 



terized in terms of the Fermi wave number, p = p 3 F /3ir 2 . The gap at zero temperature and 
the critical temperature for unpairing are listed in the Table HI The dependence of ratio 
A(T)/T on T/T c in this model is non-universal, i. e. it depends on the density; this is 
contrary to the prediction of the BCS theory with contact pairing interaction. 

Within the adopted Fermi-liquid description of neutron matter, the particle-hole interac- 
tion is given by 

VP h (q) = v v (q) + v A (q)(*.*'), (14) 

where cr refers to the Pauli matrix. The Landau parameters v v (q) and v A (q) depend on 
the momentum transfer q in the process where both fermion momenta are on the Fermi- 



surface 35( . The tensor and spin-orbit terms are small in neutron matter and were neglected 
in Eq. ffT4j) . The full dependence of these parameters on the momentum transfer is commonly 
approximated by a Legendre polynomial with respect to the angle formed by the incoming 
fermions, whereby only the leading and next-to-leading order terms contribute significantly. 
Table [I] lists the effective mass, the zeroth order Landau parameters in the scalar and 



spin channels computed within the formalism of Ref. 36J starting from the CD Bonn poten 



tial |37(]. The solution of the gap equation was obtained by applying the iterative method 
with "running" cut-off 1 38] whereby the effective pairing interaction was approximated by 
the Gogny DS1 force |39l |. 
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N N 
FIG. 2: The one-loop contribution to the polarization tensor in the superfluid matter; solid 

lines refer to the baryon propagators, wavy lines to the (amputated) Z° propagator. 
III. THE PAIR-BREAKING PROCESSES AT ONE-LOOP 

The neutrino emissivity (the power of the energy radiated per unit volume in neutrino- 
anti-neutrino pairs) is given by j^, Q] 

(G \ 2 f d?q\ f d s q 2 f 
U2j J Wyk^J j2^k-J dqo 

x Jd 3 q5(q 1 + q 2 -q)5(uj 1 + uj2-qo)qo 

x g(q )A^(q u q 2 )^mU K (q), (15) 

where G is the weak coupling constant, qi = i — 1,2 are the on-mass-shell four- 

momenta of neutrinos, g(qo) = [exp(q /T) — is the Bose distribution function, ITx(g) is 
the retarded polarization tensor, A MA (g 1; q 2 ) = Tr [7 M (1 — 7 5 ) j/i7 w (l — 7 5 ) i 2 \ Here and be- 
low the emissivities are given per neutrino flavor; the total rate is obtained upon multiplying 
the single flavor rate by the number of neutrino flavors within the standard model, Nf = 3. 
The polarization tensor of baryons at one-loop is shown in Fig. [2] and is given analytically 
by 

Il V l A (q) = TJ2 [G( P )G(p + q) T F(p)F*(p + q)] , (16) 

where p = (ipo,p). The upper/lower signs correspond to vector current (V) and axial vector 
current (A) couplings. In writing Eq. ffToT) we assumed that baryons carry the same isospin 
quantum number. 

Performing the Matsubara sums in Eq. (TT6]) (see the Appendix) we obtain 

= J2 [/(<*) " /(<*)] (~-^ ■ \ ) 

^ \iq + e p -e k iq-e p + e k J 

+ E ^ H*) " ■ ,\ ) , (17) 

\iq-e k -e p iq + e p + e k J 

where k = p + q, A T = u 2 p u\ =F h, B T = v%o\ =F h, C T = u\v 2 v ± h, D T = u 2 v v\ ± h, 
h = u p u k v p v k . The emissivity ffToT) requires the imaginary part of the polarization tensor 

8 



which after analytical continuation in Eq. (ITT)) becomes [note that /(— e p ) = 1 — f(e p )] 



Smn^(?) = - n J2if(tp)-f(tk)](A T + B T ) 6{u + e p -e k ) 

up 

- ttJ2 lf(-£ P ) - /(<*)] % " ep - <*) " % + e P + e fc )] . (18) 



To obtain the first line we used the fact that the quasiparticle spectra are invariant under 
spatial reflections, i. e., e(—p) = e{p). 

The first line in Eq. (ITS)) corresponds to the process of scattering where a quasiparticle is 
promoted out of the condensate into an excited state, or inversely, an excitation merges with 
the condensate. The corresponding piece of the response function QmH v ^ A (q) vanishes for 
small momentum transfers. Indeed neutrino energies are of the order of temperature which 
implies that their wave vectors g[fm _1 ] ~ uj u /hc ~ T/hc ~ 1/197.3 <C 1. On the other 
hand the neutron wave vectors p ~ pp ~ 1 fm _1 . Upon expanding the argument of the 
delta- function with respect to \q\ around the point q = one finds 



uo — q 



dq 



q 



A p dA p+q 



0. 



g=0 



(19) 



9 =o " e p dq 

If we assume that A ^ A(p) the third term on the 1. h. side vanishes. It follows that 
x = (p ■ q)/pq = (e p /£p)(u/vq) < 1, where v (~ Vp) is the baryon (Fermi) velocity. For on 
mass-shell neutrinos u = cq and the latter condition can not be satisfied, i. e. the scattering 
contribution can be neglected. The non-locality in the momentum or frequency domains of 
the gap function will alter this conclusion, but will require a specific model of the momentum 
and frequency dependence of the gap function (for 5-wave interactions the momentum and 



energy dependences are described respectively in Refs. [30|, [39|, |4l|, |42j and [25|, [271, |29j). 
The second line in Eq. ( ITS)) describes the process of pair-breaking and recombination, i. e., 
excitation of pairs of quasiparticles out of the condensate, and inversely, restoration of a pair 
within the condensate. Since we are interested in the emission process we shall keep only 
the terms that do not vanish for to > 0; then, the pair-braking contribution is given by 



Smn^(g) = -^[/(-g - f(e k )] C T . 



(20) 



up 
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r 2 = ^C^.^^^* +^^^^^>^><i ^ ^>>^ 



FIG. 3: Coupled integral equations for the effective weak vertices in superfluid baryonic matter. 
The "normal" T\ vertex (full triangle) and two "anomalous" vertices T2 (hatched) and 1^ (shaded 
triangle) are shown explicitly, the fourth vertex (empty triangle) is obtained by interchanging the 
particle and hole lines in the first line. The anomalous vertices vanish in the normal state. 

In the limit q — > and assuming that A 7^ A(p) the integrations in Eq. ( 1201 can be 
performed analytically. One finds 

-SmU v (q) = -2w{p F )g{u>y l f {^f (^Q ^2! AA 2 6 ^ ~ 2A )> ^ 



Smh^g) ~ + O ( , (22) 



where p{]Pf) — m*pF/2ir 2 is the density of states (h = 1) and 9 is the Heaviside step 
function; the explicit form of the 0(v 2 F /c 2 ) contribution to the axial current response is 
given in Ref. jlj. Note the threshold behavior of the vector current response, which is finite 
for frequencies that are large compared to 2A - the energy needed to break a pair. 

and carrying out the phase-space integrals we 



Upon substituting Eq. ( I2T1) in Eq. ([15 
obtain the emissivity per neutrino flavor 



x 5 „ fx 



where z = A (T) / T and 



^ = «*L*w=w*™ (23) 



* = (24) 
60 7i 6 



IV. WEAK INTERACTION VERTICES 



Consider the weak vector and axial vector vertices in the nuclear medium featuring a 
condensate. Because the particle-hole interactions in the medium (which are represented 

10 



by the Landau parameters in Table H]) are not small, the vertex renormalization requires 
a summation of an infinite number of particle-hole loops. There are four topologically 
non-equivalent vertices in the superfluid state in general, but, as we shall see, the particle- 
hole symmetry reduces their number by one. Since neutrons pair in an isospin-1 state 
(neutron-proton pairing is unimportant for large asymmetries 43j) we shall suppress the 
isospin indices. The expansion of the Landau parameters in Legendre polynomials will be 
truncated at the leading order (the next-to-leading order terms are suppressed by powers of 
vp/c). Thus, the effective particle-hole interaction is given by 



v v + v a (<t-<t'). 



(25) 



The integral equation defining the effective weak vertices, which we write in an operator 
form, are given by 



1 1 



1 2 



1 3 



+ v a (GT a 1 G + FYIG + GV a 2 F + FT^F), 
v a (GT a 2 G ] + FTIG ] + GTIF + FT^F), 
v a (G^G + FT^G + G^TIF + FT%F), 



(26) 
(27) 
(28) 



and are displayed diagrammatically in Fig. [31 Here F = —ia y F, v a with a G V, A are 
defined in Eq. (1251) . while Tq = 1 and Tq = <r. The fourth integral equation for the vertex 
follows upon interchanging in Eq. (|26l) particle and hole propagators G <-> G'. The 
momentum space representation of Eq. (1261) reads 

d 4 p 



r?(g) = T a + v a 



{2k)< 



G(p)Y a MG(p + q) + G(p)T a 2 (q)F(p + q) 



+ F{p)Tl{q)G{p + q) + F{p)T a 4 (q)F{p + q) 



(29) 



with similar expressions for Eqs. ( 1271) and ( |28l) . Even though the driving interactions are 
local in time, the resummations at finite temperatures lead to time-retarded interactions, 
which imply that the effective vertices in Eqs. (|26l)-(l28l) are complex in general. Considering 



the scalar interaction Tq we obtain 



( (v v )- 1 



- [U GG (q) - U FF (q)] 
-2Tl GF (q) ( /■ 

-2U FG (q) 



n GF (g) 



n FF (g) 

Gte(?) 



; y )- x -n 



r 2 v (?) 



1 





V u J 



(30) 
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where, using the abbreviations X, X' e G,F,G^, we defined (see the Appendix) 



n 



XX 



(?) 



/ w? Lxx ' {q ' p) = I Wf^ x{p)x ' {p + q) - 



d 3 p 



J G\GU Lqf 



(31) 



and 



The loop integrals posses the following symmetries: L GG - 
L FG — L G tp, which imply that = Y\ and = —Tf. For energy-momentum inde- 
pendent interactions the integral equation ( 1301) reduces to three coupled linear equation for 



the complex functions T 



v 



1, 2, 3. The details of the computation of the coefficients in 



(32) 



Eqs. f l30|) are relegated to the Appendix. 

In the case of axial-vector interactions the effective vertices are given by 

U GF (q) n FG (g) 
(tr 4 )- 1 - n GGt (<?) U FF (q) 

o u FF (q) Kr 1 -n GfG ( g ) J 

This equation has a "trivial" solution T± = 1 and T 
v F /c 2 expansion and we shall not consider it further. In the weak-coupling BCS limit there 
exist approximate constraints among the loop integrals 



( (v A ) 



\ 



^ i rf( q ) ^ 

r£(g) 

V r 3 (?) J 



f T A\ 
1 



V / 



Tf = to the leading order in the 



J GF 



— L 



FG, 



J G J <G 



J GG J < i 



(33) 



which allow us to reduce the number of equations in the set (I30p from three to two. The 
relations (13"31 are exact at the threshold u = 2A and hold approximately for systems with 
strong degeneracy (corrections being suppressed by powers of the ratio of the temperature 
over the Fermi-energy). It follows then from Eq. (|30|) that T\ = — Tg . The solutions for the 
remaining vertices Ti and T2 are convenient to express through linear combinations of the 
polarization tensors (1311) . defined as 

d 3 p 



/dp 

d 3 p C p (u,q) 



(27T)* 



B(q) = 2U FG (q) = -luA 

C(q) = n GGt (q) + U FF (q)-(v v )- 1 
where (k = p + q) 



d 3 p 
(2^p 



[2e p L p (0,0)-F+(q)£ p (uj,q 



£ p (uJ,q) 



L£k A 2 



f{€p)-f(*k) 

(e p + e k ) 2 + iS' 



(34) 
(35) 
(36) 

(37) 
(38) 
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n( a ) 

FIG. 4: The sum of polarization tensors that contribute to the neutrino emission rate. The 
contributions form H^(q) and U^(q) vanish at one-loop. 

In obtaining Eq. fl36l) we used the fact that 

1 + vV / A T^^GGt + ILff) = 0, (39) 

J y^ 71 ) w=0,q=0 

where A is a three-dimensional ultraviolet cut-off on the momentum integration, which is 
required for regularization of the gap equation fl39l . A may be adjusted to reproduce the 
gaps obtained from finite range interactions in Sec. [Til but such an adjustment is not required 
as long as the constraint ( |39i) is taken into account in finding the solutions for the vertex 
functions. 

The vertex functions expressed through the combination ( !34j) -( l36j) read 

Tliq) = C(q)-vV[A(q)C(q) + B(q)*]' m 

r / x B{g) , , 

2{q> ~ C(q)-vV[A(q)C(q) + B(qyY 1 ' 

The poles of the functions Ti(q) and ^(g) determine the collective excitations of the S wave 
superfluid, which we will discuss in the next section. 



V. THE FULL POLARIZATION TENSOR AND THE EMISSIVITY 

Having determined the effective vertices T^o;) and T^(u;), we now construct the com- 
plete polarization tensor, which is given by the sum of diagrams in Fig. HI Summing the 
contributions we obtain the full polarization tensor expressed through the combinations 
(El- (EE}: 

U V (q) = A{ q )C{q) + B{qf 

{q) C(q)-vV[A(q)C(q) + B(qyY 1 ] 

Eq. f|42|) is our central result valid for arbitrary momentum transfers. It can be used to 
compute the emissivity directly, but it is illuminating to work with its small q expansion. 
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The leading order term in this expansion vanishes, i. e. the full polarization tensor vanishes 
when q = (note that this is contrary to the one-loop polarization tensor, which is finite 
in this limit, see Eq. (1171) and the following discussion.) Indeed, taking the limit q — > in 
Eqs. AMD- flMD and substituting in Eq. (02} we find that 

]imU v (q,u) = 0. (43) 



This result is consistent with the /-sum rule 44j 



lim / duj uj lmll v (q,uj) = 0, (44) 

which follows directly from (|43p . The opposite need not to be true for arbitrary functional 
dependence of ImU v (q, uj) on frequency, but is the case in practice. The reason is that the 
causal polarization tensors are odd functions of frequency. Furthermore, since they should 
correspond to stable collective modes the condition uj \mH v (q, uj) > is satisfied, which 
combined with the /-sum rule leads back to Eq. (j43|) . 

Consider now the next-to-leading order terms. Keeping only the leading order ^^(0) 
function in Eqs. ( I34l) - (l36l) one finds that the polarization function vanishes order by order in 
the expansion of the function £ p (uj, q); thus instead of straightforward expansions of kernels 
in Eqs. f l34|) - fl36|) . we shall expand only the functions J-'p(q): 

J?(q) = ( ejL p) (l + ^ + ^-)=2e p + ^ (45) 
p V 2 I V e p e k e p e k J e p 



T - M = to) U - &k + ^1) = ^ - ( 46) 



and 



tp^k £p^k / £p £p 

where £ q = q 2 /2m is the nucleon recoil (the linear in q terms are omitted since they vanish 
upon angle integrations). Substituting these expressions back into Eqs. ( !34l) - (|36l) one finds 

A(q) = 2A 2 I (q)-A%I A (q), (47) 
B{q) = -ujA I (q), (48) 
C(q) = -y I (q)+^Ic{q), (49) 
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where we defined the following integrals 



d^p £ 



(2tt) 3 e 



(50) 
(51) 
(52) 



The pole(s) of the polarization tensor (or equivalently of the vertex function) determine the 
dispersion relation of the collective excitations: 



V{q) = C(q)-v v [A(q)C(q) + B(q) 2 } 

~ Co(q) + &(q) - v v [A 1 (q)C (q) + Mv)Ci(q)] = 0, 



(53) 



where the subscripts and 1 refer to the leading and next to leading order terms in the small 
q expansion. Substituting here the expressions (H7I) - (HU1) we obtain the dispersion relation 
for the acoustic modes 

u 2 = c 2 q 2 , (54) 

where the sound velocity is defined as 

Ic 



2m* 



r-1 



4v v A 2 



l+UJ 



2 I A 



4L 



c 



(55) 



When v v = 0, the sound velocity at zero temperature should reduce to that of non- 
interacting Fermi-gas, c 2 = Vp/3. It is seen that the ratio Ic/Iq ~ £ P ; since the average of 
£p ~ VfVf one recovers the scaling c 2 ~ v 2 F . 

At the next-to-leading order in small-g expansion the polarization tensor is given by 



UX{q)=Ai(q) + 



Co(q) 



cm = 



I A (q) + —I C {q) 



(56) 



It is sufficient to evaluate the integrals in the limit q = at the order we are working; then 
the imaginary part of the polarization tensor can be evaluated analytically 



ImlT^) = -8nu(p F )-^g(iu)- l f (|) 6(lu - 2A). 
The neutrino emissivity to the leading 0(q 2 ) order is 

e^p = e z 2 ^- I dx x 5 f ( - 



(57) 



(58) 
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FIG. 5: (Color online) Dependence of neutrino emissivity in units of eo for one-loop (upper panel) 
and full (lower panel) polarization tensor on the reduced temperature for Fermi wave vectors 
Pi?[fm~ x ] = 0.4 (solid), 0.8 (dashed, red online), and 1.4 (dashed-dotted, green online). 

Compared to the one-loop result, the emissivity is suppressed by a factor T/m* and shows a 
different functional dependence on the frequency; this difference arises from the parametric 
suppression of the full polarization tensor IT] 7 (q) ~ £ 9 and the fact that in the bremsstrahlung 
process q ~ T. The emissivity is thus suppressed compared to the one- loop result, roughly 
by a factor 4 x 10 -3 , as illustrated in Fig. [5j Note that the density dependence in Fig. [5] 
arises from the density dependence of the ratio A(T)/T as a function of T/T c (in the BCS 
theory with contact interaction the ratio e uP /eQ is universal). For T — > T c the both rates 
vanish, consistent with the observation that the pair bremsstrahlung is absent in normal 
matter for on-shell baryons. At small T < 0.3T C the rates are suppressed exponentially as 
exp(— A/T). At intermediate temperatures the emissivities that include vertex corrections 
should be scaled down from their one-loop counterparts roughly by the factor ~4x 10 -3 
quoted above. 
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VI. SUMMARY 



The neutrino losses via neutrino pair bremsstrahlung from neutron superfluid at subnu- 
clear densities have been studied within a microscopic model based on the Fermi-liquid and 
BCS theories. In our model the modifications of the weak interaction rates in the nuclear 
medium are taken into account by summing infinite series of irreduceable particle-hole dia- 
grams in terms of a contact (momentum and energy independent) driving interaction. These 
modification are embodied in the vertex functions that are computed from polarization ten- 
sor components corresponding to the pair-breaking process within the kinematical domain 
uj G [2A; oo] and small-g expansion. Special attention was paid to preserving the dispersion 
relations for the polarization tensor components, that are implied by the unitarity of the 
^-matrix for the emission process. The leading order contribution to the rate arises from 
0(q 2 ) contribution to the polarization tensor, which vanishes identically when q = 0. We 
find neutrino loss rates that are suppressed compared to the one-loop results by a factor 
of the order 4 x 1CT 3 . The magnitude of the suppression differs from the one predicted in 
Ref. |l5| due to two factors: (i) the vertex functions used herein include the finite temper- 
ature effects, which guarantees their consistency with the associated polarization tensors; 
(ii) the full polarization tensor is derived in a different expansion (small-g rather than the 
v F jc expansion). 

The modifications to the neutrino emission rate through pair-breaking process found 
above call for a detail reassessment of their role in the late-time cooling of neutron stars. 
Even though the rates are suppressed substantially, their strong dependence on the value of 
the pairing gap does not rule out the possibility that for large enough gaps the rates will 
become comparable to those of the competing processes. 

While we concentrated above on the neutral current interactions, our formalism can 
be adapted to compute the rates of charge-current Urea process n — * p + e + v e beyond 
the one- loop rate 45]. These should be subject to strong correction due to the tensor 
neutron-proton correlations. Vertex corrections could be important in the related problem 
of neutrino emission and propagation in quark matter 14j| , where one-loop results for a 



number of processes became available recently [46|, |47|, 



48, 



49 



50|,|5l|. 
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APPENDIX A: LOOP INTEGRALS 



Here we quote the results for the loop integrals and polarization tensors that have been 
used in the main text. The loop integrals are defined as convolution products of the Mat- 
subara Green's functions (here k = p + q) 



L GG {q,P) = Tj2G{tp,p)G{tp + tq,k) 



2 2 



2 2 
V p V k 



iq + e p - e k iq - e p + e k 



+ 



2 2 
U k V p 



2 2 
U p V k 



iq-e p -e k tq + e p + e k 
L FG (q,P) = Tj2 F ( l P,P)G(ip + iq,k) 



lf(e P ) - /(£*)] 



lf(-e p ) - f(e k )} , 



Ip 



UpVp 



- + -. — K - — 

iq + e p - e k %q - e p + e k 



+ u p v p 



ui 



+ 



iq -e p - e k iq + e p + e k 
L FF (<1,P) = TY,F(ip,p)F\ip + iq,k) 



lf(-e P ) - f{e k )\ , 



u p u k v p v k 



+ 



iq + e p - e k iq - e p + e k 
1 



lf(e P ) ~ f(e k )} 



iq + e p + e k iq - e p - e k 



[f(-e P )-f(e k )]}, 



L& G (q,P) = TY,G\ip,p)G(ip + iq,k) 



%p 



u 2 p v 2 k 



iq + e p - e k iq - e p + e k 



u 2 p u\ 



VpVl 



iq- e p - e k iq + e p + e k 



Wp) - f(e k )} 
lf(-e P ) ~ f(e k )} , 



(Al) 



(A2) 



(A3) 



(A4) 
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L FG t{q,P) = T^F{i V ,p)G\ip + iq,k) 



UpVp 



UpVp 



(II 



-. — K - — + - 

tq + e p - e k iq - e p + e k 



+ 



in 



[/(ep) - /fa)] 
[f(-e p ) - f(e k )} . 



(A5) 



_iq-Ep-e k tq + e p + e k 
The remainder loop integrals are obtained from those defined above through the relations 



L&Gt(iq,p) = L GG (-iq,p), L GF (iq,p) = L FG (-iq,p), 
L GG i(iq,p) = L G i G (-iq,p), L GtF (iq,p) = L FG t(-iq,p), 



(A6) 
(A7) 



where, except for the first relation, we used the fact that the quasiparticle spectrum is 
reflection invariant, e(— p) = e(p). This property implies that the arguments of the functions 
can be interchanged k «-> p without changing the result. Furthermore, upon performing the 
substitution p — > — p — q one finds that 

L GG (iq,p) = L G f G f(iq,p), L GF {iq } p) = L FGt (iq,p), L FG (iq,p) = L GtF (iq,p). (A8) 

The retarded polarization tensor is obtained by analytical continuation in the loop integrals 
Lxx' (iq, p) = L X x' + id, p) and by subsequent integration over the three- momentum p 
according to the Eq. (!3T|) . 
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